Persistent currents in Mobius strips 
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Relation between the geometry of a two-dimensional sample and its equilibrium physical properties 
is exemplified here for a system of non-interacting electrons on a Mobius strip. Dispersion relation 
for a clean sample is derived and its persistent current under moderate disorder is elucidated, 
using statistical analysis pertinent to a single sample experiment. The flux periodicity is found to 
be distinct from that in a cylindrical sample, and the essential role of disorder in the ability to 
experimentally identify a Mobius strip is pointed out. 
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I. INTRODUCTION 

An important theme in quantum mechanics is to find 
a relation between a global geometry of a sample (e.g., 
boundary conditions) and its physical properties. We 
address this issue by comparing flux periodicity of per- 
sistent currents in a cylinder and in a Mobius strip. The 
aim is to determine whether the geometrical (in some 
sense topological) difference is tangible and experimen- 
tally observable. At zero temperature, the persistent cur- 
rent I ((f)) in a ring can be expressed asi 



dE((f>) 



I n sin(27m</ 



(1) 



where (f> is the magnetic flux threading the ring in units 
of $o = hc/e, E(<f>) is the ground state energy, and J„ 
are the current harmonics. 

The current I ((f)) is an anti-symmetric and periodic 
function of <f> with period 1. Possible occurrence of 
smaller flux periodicity in mesoscopic physics is one of the 
cornerstones of weak localization. For the cylinder geom- 
etry, conductance measurements^ and magnetization of 
10 7 copper rings^ indicate the emergence of periodicity 
1/2. It is shown to be intimately related to the procedure 
of averaging over disorder realizations and numbers of 
electrons in the ringsS^^. Very recently, a microscopic 
NbSe3 Mobius strip has been fabricated 7 -. Obviously, in 
this case, attention should be focused on a single sample 
measurement^ for which there is no self-averaging. 

Let us first mention several intuitive points relevant to 
the flux periodicity in the Mobius strip, based on semi- 
classical arguments and geometry 9 . First, recall that the 
periodicity is related to interference between trajecto- 
ries (such as Aharonov-Bohm interference between differ- 
ent trajectories or weak-localization interference between 
time-reversed paths). In the cylinder (Mobius) geometry, 
an electron moving in the longitudinal direction along the 
ring encircles the system once (twice) before returning to 
its initial position. Therefore, we might expect differ- 
ent flux periodicities of the persistent current between 
the two cases. Second, unlike a cylinder which can be 
"pressed" into a one-dimensional ring, the Mobius strip 
cannot be pressed into a one-dimensional structure. This 
brings in another important factor, namely, the motion of 
electrons in the transverse direction. In a tight-binding 



model this motion is controlled by the transverse hop- 
ping. If it is very weak, the twice-encircling property of 
the Mobius strip implies the dominance of even harmon- 
ics I 2n - Contrary, for a strong transverse hopping, the 
current in the Mobius strip is expected to be effectively 
similar to that in the cylindrical strip 10 . In the following 
we are mainly interested in a regime where the trans- 
verse hopping is slightly less than with the longitudinal 
one. Third, the role of disorder should be carefully exam- 
ined. Weak disorder is not expected to significantly alter 
interference between semi-classical trajectories discussed 
above, while strong disorder should result in a reduced 
sensitivity to the pertinent geometry, due to localization 
effects. The most intriguing disorder effect might then be 
expected in a moderate strength of disorder which will 
be used below. The upshot of the present study is that 
the periodicity pattern in a Mobius strip is remarkably 
distinct from that of a cylinder, and that disorder plays 
a crucial role in making the statistical effect detectable. 



II. MODEL 

A Mobius strip is modelled by considering a non- 
interacting particle in a rectangle of length L x and width 
L y , requiring its wave-function if>(x,y) to satisfy Dirich- 
let boundary conditions in the y direction, and Mobius 
boundary conditions^ in the x direction: 



^(x,~L y /2) 
ip(x+L x ,y) = 



= V(z, L y /2) = 
ip(x,-y) 



(Dirichlet B.C.) , (2) 
(Mobius B.C.) . (3) 



The quantized wave-numbers are k v = (Tr/L y )n y 



n x ), where 



1,2, 



and 
and 



- (2vr/^)([| 

n x — 0, ±1,±2, The notation [a) n represents a for 
n = even and for n — odd. In the cylinder geometry, 
Eq. (3) should be replaced by yj(x + L x , y) = if)(x, y), and 
gives k x = (2tt/ L x )n x . Thus, only the n y = even eigen- 
states are affected by the switch from the conventional 
cylinder (periodic) boundary conditions to the Mobius 
ones. 

In the absence of disorder, the energies of the eigen- 
states both in the Mobius and cylinder strips are given 
by the formula 



p. 



k x -^)+e y (k y ), (4) 
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where e x and e y provide the dispersion relation. Equa- 
tion is rather general for clean systems. To be more 
specific, let us model the Mobius strip by a tight-binding 
Hamiltonian. The Mobius strip is constructed from a 
rectangular lattice including N x 2M sites. The rectan- 
gle is twisted by 180°, and its two sides are connected, 
such that longitudinal wire 1 is attached to wire 2M, 
wire 2 is attached to wire 2M — 1 and so on. The Mobius 
strip so constructed includes M longitudinal wires with 
2N sites on each one. The Hamiltonian is then 



2N M 



H 



Mobius — ^ ^ ^ ^ ^rim^mCnm t\C ^ C^ m C 7? ,-|-i r 
n— 1 in — I 
2N M-l 

C nm+l C n 



2i\ 



c \ lM C n+NM + h.C. (5) 



where c nm is the fermion operator at the site (n, m) 
(n = 1,2,..., 27V, m = 1, 2, . . . , M) and h and t 2 are 
longitudinal and transverse hopping amplitudes respec- 
tively. The quantity e nm is the site energy. Connecting 
the two sides of the rectangle without twisting, we obtain 
a cylindrical strip which includes 2M longitudinal wires 
composed of N sites. The Hamiltonian of the cylinder is 



N 2M 



H, 



cylinder 



= Wfe c f c -Ue' 2m ' t ' /N c i 

— / J / J \c-nm<- nm (-nm '1 B ^ n 
n—1 m—1 

N 2A/-1 

■h^2 C nm+l C «™ + h ' C - 



nm C n+lm 



n—1 m—1 



Locally the two Hamiltonians JSJ and ijfjj look the same. 
But there is a couple of essential differences between 
them: a) The Mobius Hamiltonian JSJ includes an extra 
term which describes long range hopping between distant 
parts of the Mth wireii. b) While the magnetic phase 
accumulated along the longitudinal direction on each link 
is the same (that is, 2w<p/N), the corresponding number 
of links is different (2N for the Mobius strip and N for 
the cylinder). 



III. THE SPECTRUM 

We first consider a system without disorder, namely, 
£nm = 0. The dispersion relation for an electron in the 
Mobius strip reads, 



2tt f 


T 




.2. 



2t\ cos 



2t 2 cos n, 

\2M+1 l 



(J) 



where n x = 1, • • • , N and n v = 1, • • • , 2M . Defining new 
indexes k = [l] n + 2n x and q — [i]^ + n y /2, one obtains 
a more suggestive form, 



E, 



2t\ cos 
2^2 cos 



N 



[k - 2(f) 



2MTl {2q - [l]k) 



(8) 
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FIG. 1: Single-particle energy spectrum as a function of flux 
threading the Mobius ring. The 195th-200th energy levels are 
shown. The parameters are = 20, M — 10, = 0.5, and 
W = 0.5. Energies are measured in units of ti. 
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where k = 1, • • • , 2N and q = 1, • • • , M . It is instruc- 
tive to compare it with the energy in the cylinder geome- 
try, E kq = ~2t x cos [22 (k -(/))] -2i 2 cos (jM+il)' where 
k = 1, • • • , N and q = 1, • • • , 2M. Despite the apparent 
similarity between these two spectra, there are at least 
two important differences. First, the combination of flux 
and longitudinal momentum is distinct, namely, it is k— <p 
for the cylinder and k — 2(f) for the Mobius strip. For a 
small ratio ^Ai this might affect the periodicity of the 
currenljiS. Second, the mini-band structure is different. 

We now turn to elucidate the current in disordered 
Mobius strips. The random numbers e nm are assumed to 
be uniformly distributed over the range —W/2 < e nrn < 
W/2, where W represents the strength of disorder. The 
Hamiltonian Eq. JSJ) [or Eq. ©] is treated numerically. 
As an example, the evolution of single-particle energies 
with flux in a disordered Mobius strip with N = 20 
and M = 10 is shown in Fig. 1. The parameters are 
t2/ti = 0.5 and W/ti — 0.5. The pattern of avoided 
crossing turns out to be remarkably different from that 
for a cylinder (see ref. 5 figure 1 therein). It must then 
be reflected in the behavior of persistent currents. 

The first stage of the analysis is an inspection of the 
typical values of /„, aiming in determination of their de- 
pendence on the ratio tijt\. As expected, in the absence 
of averaging we find typical 1\ dominance in case of the 
cylinder geometry irrespective of the tijt\ ratio. For the 
Mobius geometry the emerging picture is quite different. 
Figure 2 shows the Fourier components of the persistent 
current for a clean Mobius strip as a function of the ratio 
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FIG. 2: Fourier components of the persistent current for the 
clean Mobius strip as a function of the transverse-hopping 
energy (a) at the half-filling and (b) below the half-filling. 
The size of the Mobius strip is given by N = 20 and M — 10. 



t 2 /ti at and below the half-filling (N e — 200). For small 
ratios (ti/tx < 0.1) we find, as can be naively expected, 
I2 dominance. The expected effect of averaging in the 
cylinder case is to emphasize the 1% contribution, while 
in the Mobius case the expected effect is to emphasize 
the I4 contribution. For clean Mobius strip the /„ with 
odd n, as a function of the number of electrons N e , is 
anti-symmetric around half-filling. Therefore I\ and 1% 
completely vanish [Fig. 2(a)]. See further discussionin 
Sect. V. To avoid this particularity at the half-filling, we 
display in Fig. 2(b) also the case where the number of 
electrons (N e = 190) is below half-filling. For large ra- 
tios {t 2 /t\ > 0.8) we observe in Fig. 2(b) a cylinder-like 
regime where there is typically I\ dominance. This is be- 
cause the strong transverse hopping changes the period- 
icity of the Mobius strip to that of the conventional cylin- 
der. The somewhat unexpected observation is that there 
is a distinct wide intermediate regime (0.1 < tz/tx < 0.8) 
where I\, I2, I3 and I4 are all comparable. This is the 
regime which is of experimental relevance. The expected 
effect of averaging in this regime is to emphasize both 
the I2 and the I4 contributions. 



IV. STATISTICAL ANALYSIS 

The problem arising in the analysis of persistent cur- 
rents in disordered Mobius strips is how to character- 
ize the statistics of the calculated data. It was already 
pointed out that essential properties of observables result 



from the averaging procedure and the nature of the un- 
derlying statistical ensembl o 4 i 5 i 6 . On the other hand, fab- 
rication of a Mobius strip requires an outstanding effort 7 , 
and hence, anticipated measurements of the persistent 
current would probably be performed on a single sam- 
ple. Thus, somewhat unfortunately, the important re- 
sults reported therein and the powerful calculation meth- 
ods based on super-symmetry might be less useful for 
single-sample experiments since there is no averaging. 

What is then the most efficient way to present our 
calculated results? The answer is provided by elemen- 
tary statistics. An experimental result consists of a set 
of K measurements I((f>i), i = 1,2, ... , K performed on 
a given sample. This sample is taken out of an en- 
semble of Mobius strips with different disorder realiza- 
tions, electron numbers N e , aspect ratios, etc. The set 
{I((f)i)\i = 1, . . . , K} can be regarded as an instance of a 
random vector in a it - dimensional space. Alternatively, 
this instance can be represented by the current harmonics 
(7i, I2, . . . ) defined via Eq. QJ. For our purpose it seems 
adequate to keep only the first 4 harmonics. The relevant 
statistical ensemble is then a set of "points" (I\, I2, 13,14) 
in four-dimensional probability space, each point corre- 
sponds to a possible experimental measurement of the 
current on the entire (f> interval. Let us denote the num- 
ber of points within an infinitesimal four-dimensional vol- 
ume element by P{I\, I2, I3, I^dlid^dl^dl^. The distri- 
bution function P is normalized to Af, the total number 
of members in the ensemble. The most probable (typical) 
experimental result is then determined by the quadruple 
Iii 12,13, 1 a at which P is maximal. Another quantity, 
which seems more informative and easy to analyze, is 
the distribution 



Pn(In)= P{h,h,h,h) J] d\I m \. (9) 



This corresponds to the possibility of finding a sample 
whose current /(</>) is approximately described by /(</>) ~ 
I n sin(27m0). (For a sample counted by p n (I n ), all the 
harmonics other than I n are at most half of in mag- 
nitude). The number of members in the ensemble that 
exhibit /„ dominance is therefore Af n = p n (I n )d\I n \. 
If Af n > AT m for any m ^ n, the typical periodicity of 
I(4>) is dominantly 1/n. In actual calculations, we as- 
sume that the lattice structure, the aspect ratio, and the 
strength of disorder are fixed, and that the temperature 
is very low. Then, two quantities are still fluctuating, 
namely, the filling factor (or the electron number N e ) 
and the specific realization of disorder. We generate an 
ensemble of Af = Af a Af b members corresponding to Af a 
consecutive values of N e , usually around half filling, and 
Af b realizations of disorder for each one of them. Actu- 
ally, for our systems of size N = 20, M = 10 with t\ = 1, 
t 2 = 0.5, and W = 0.5, we take 150 < N e < 250, hence 
N a = 101 and N b = 250, so that N = 25250. The dis- 
tributions Pn(I n ) for the cylinder and Mobius ensembles 
are shown in Fig. 3. 
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FIG. 3: The distributions p„(I n ) (n = 1,2,3,4) defined by 
Eq. JjJJ for the cylinder and Mobius ensembles. The numbers 
of members with dominance are A/*i = 15, 829, A2 = 382, 
A3 = 4, and A/4 = 439 for the cylinder ensemble, and A/*i = 
1, 562, A/2 = 336, A/3 = 384, and Ha = 1, 992 for the Mobius 
ensemble. 



FIG. 4: 7i and I2 as a function of N e for the ordered (solid 
Une) and disordered (dots) systems. Parameters describing 
the systems are the same with those for Fig. 3. 



V. MAIN OBSERVATIONS 

The most striking result that can be deduced from 
Fig. 3 is the essential reduction of A/i for the Mobius en- 
semble compared with the cylinder one. For the present 
ratio N/2M = 1, there is also a strong tendency towards 
$o/4 periodicity, since A/4 > Af m ^4 for the Mobius en- 
semble. This result is intriguing, because here we have 
no averaging procedure which is crucial to get the 1/2 
periodicity in cylindrical strips. However, this 1/4 peri- 
odicity emerges only for the specific ratio N/2M — 1. We 
have calculated the distributions p n (I n ) for Mobius strips 
with several aspect ratios. The value of N n depends on 
the aspect ratio. No specific n gives prominent Af n in- 
dependently of the aspect ratio. On the other hand, the 
collapse of 1% dominance in the Mobius ensemble is ro- 
bust and persists in systems with different ratios N/2M 
as well. We can safely say that A/i, A/2, A/3, and A/4 
become all comparable in the Mobius ensemble. 

The natural question that comes to mind is whether 
this result is a consequence of the Mobius geometry or, 
rather, it is due to the presence of disorder. In order 
to answer this question, we have performed the calcula- 
tion of P°(Ii, I2, 13, 14) for a "clean" Mobius ensemble 
(without disorder, only N e is being changed). We found 
out that the probability to find any /„ dominance is ex- 
tremely small. The immediate conclusion is that disor- 
der is essential for the identification of Mobius strips via 
I n >i dominance. Does this mean that interference or 



weak-localization effects due to the presence of disorder 
is important? To clarify this point, we should understand 
how P°(Ii, I2, 13, 14) is modified by disorder. The distri- 
bution P (Ji, I2, 13, 1 a) is, in fact, a function defined on 
a one dimensional curve \I\{N e ),l2{N e ),I^{N e ),Ii{N e )\ 
in 72)^3) I4) space. For this reason, it is unlikely to 
find a sample where one of the I n is dominant. The ef- 
fect of disorder is to give some "thickness" to this curve 
(see Fig. 4). Taking into account that the amplitudes of 
In{N e ) for Mobius strips are all comparable, the thick- 
ness gives a finite probability to find samples where one 
of the is dominant. On the contrary, in the case of 
cylindrical strips, the amplitude of I\(N e ) is overwhelm- 
ingly larger than those of I n ^i(N e ), which makes it un- 
likely to find I n ^i(N e ) dominated samples even if we take 
the statistical effect of disorder into account. We should 
note here that the function /„ with odd n for the clean 
Mobius strip is an even function around the half-filling 
(N e = 200) and an odd function for odd n, while the 
function I n for arbitrary n is an even function in the 
cylinder case. 

Our findings regarding Af n for the Mobius ensemble are 
based on the fact that the amplitudes of I n (N e ) are all 
comparable for Mobius strips. As we have observed in 
Fig. 2, this is a robust statistical property in the inter- 
mediate regime 0.1 < £2/^1 < 0.8. The choice ^Ai = 0.5 
above, provides typical results for p n (I n ) and N n in case 
that t^/ti is within this distinct regime. 
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VI. CONCLUSIONS 

We have studied the persistent currents of non- 
interacting electrons in Mobius strips. The spectral prop- 
erties for a clean system were found analytically, and the 
effect of disorder on the currents was analyzed numeri- 
cally. We have found that disorder is quite essential for 
the identification of Mobius strips. The issue of disorder 
averaging is not relevant for single sample experiments, 
and hence, special care is required for statistical analysis 
of the current harmonics. The fingerprint of the Mobius 
geometry is an enhanced probability to find samples in 



which I n , with n > 1 dominates. This should be con- 
trasted with the case of cylinder geometry, where there is 
a clear ii dominance. The above assertion regarding the 
fingerprint of the Mobius geometry is correct provided 
the effect of disorder is properly taken into account. 
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